arXiv:astro-ph/0211237vl 12 Nov 2002 


Astronomy &: Astrophysics manuscript no. H3905 February 2, 2008 

(DOI: will be inserted by hand later) 


XMM-Newton observations of ULIRGs I: A Compton-thick AGN 

in IRAS19254-7245 

V. Braito^’^, A. Franceschini^, R. Della Ceca^, P . Severgnini^, L.Bassani^, M. Cappi^, G.Malaguti^, 
G.G.C. Palumbo®, M. Persic®, G. Risaliti^’® and M. Salvati^ 

^ INAF — Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy. 

^ Dipartimento di Astronomia, Universita di Padova, Vicolo dell’Osservatorio 2, 35122 Padova, Italy. 

® INAF — Osservatorio Astronomico di Brera, Via Brera 28, 20121 Milano, Italy. 

lASF — CNR, Sezione di Bologna, Via Gobetti 101, 40129 Bologna, Italy. 

® Dipartimento di Astronomia, Universita di Bologna, Via Ranzani 1, 40127 Bologna, Italy. 

® INAF — Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 34131 Trieste, Italy 
^ INAF — Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy. 

® Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA. 

ReceivediXX; Accepted: XXX 

Abstract. We present the XMM-Aeuiton observation of the merging system IRAS 19254-7245, also known as The 
Superantennae, whose southern nucleus is classified as a Seyfert 2 galaxy. The XMM-Aewton data have allowed 
us to perform a detailed X-ray imaging and spectral analysis of this system. We clearly detect, for the first time 
in this system, a strong [EW ~ 1.4 keV) Fe emission line at 6.49 ± 0.1 keV (rest-frame). The X-ray spectrum 
requires a soft thermal component {kT ^ 0.9 keV; I/ 0 . 5-2 ^ 4 x lO'^'^ erg s“^), likely associated with the starburst, 
and a hard power-law continuum above 2 keV (observed I/ 2-10 ~ 4 x 10"*^ erg s”*^ ). We confirm the flatness of 
this latter component, already noted in previous ASCA data. This flatness, together with the detection of the 
strong Fe-Ka line and other broad band indicators, suggest the presence of a Compton-thick AGN with intrinsic 
luminosity ^ 10^"* erg s“^. We show that a Gompton-thick model can perfectly reproduce the X-ray spectral 
properties of this object. 
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1. Introduction 

The Ultra Luminous Infrared Galaxies (hereafter 
ULIRGs) are an enigmatic class of sources which emit 
most of their energy in the far-IR (FIR) domain, with lu¬ 
minosities in excess to 10^^ Lq, (i.e. comparable to QSO 
luminosities). While it is now clear that the huge observed 
FIR luminosity is due to the presence of large amounts of 
dust reprocessing the “primary” optical and UV emission 
from the central source, the nature of the latter is still 
debated. It is now widely accepted that both starburst 
(SB) and AGN activity may be responsible for the ob¬ 
served luminosity, but their relative contributions are still 
unconstrained (in some cases even the presence of an AGN 
is unclear). 

Hard X-ray {E >2 keV) observations, less affected by 
the photoelectric absorption, are a fundamental tool not 
only to investigate the presence of an hidden AGN, but 
also to study its physical properties and to estimate its 
contribution to the high observed FIR emission. Indeed 
some of the ULIRGs, classified as pure SBs based on op¬ 


tical and IR spectroscopy, show spectral properties typi¬ 
cal of obscu red AGNs whe n observed in hard X-rays (i.e., 
NGC 6240: Iwasawa I99S ). 

To shed light on this topic we are carrying 
out a mini-survey with XMM-Newton of 10 nearby 
(z < 0.2) ULIRGs (see Braito et al. 2002 ) for 

which high-quality mid-IR and optical spectroscopic 


data are available (Genzel et al. 1998; Lutz et al. 1999; 
Veilleux et al. 19^ ) . XMM-Newton, thanks to its wide en¬ 
ergy band (0.2-10 keV) and high sensitivity (as compared 
to previous missions such as ASCA and ROSAT) allows 
for the first time the investigation of the X-ray spectral 
properties of this class of sources, which are usually X- 
ray faint ( Risaliti et al. 200*^ ). In this paper we discuss 
the 0.2-10 keV XMM-Newton spectrum of IRAS 19254- 
7245 ( The Superantennae), a well known ULIRG (Lfir = 
1.1 X IO^^Lq) located at redshift of z = 0.062. In par¬ 
ticular, we will argue that IRAS 19254-7245 harbors a 
high-luminosity {L ^ lO'^^ erg s“l), Compton-thick AGN. 
In this letter we assume Hq = 50 km s“^ Mpc“^ and 
qo = 0.5. 
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2. IRAS 19254-7245: Main Properties 


Like most ULIRGs studied so far, IRAS 19254-7245 
is a merger system. It shows giant tails extending out 
to a distance of 350 kpc, triggered by the merger of 
two gas-rich galaxies with nuclei ~ 9 arcsec (with 
Hq = 50 km s“^ Mpc“^ this corresponds to ~ 15 
kpc ) apart from each other. In the southern nucleus, 
which is the dominant source at different wavelengths, 
mid-IR spectrosc opy reveals the presence of an AGN 
( Lutz et al. 1999|), classified as a Seyfert 2 from o ptical 
spectroscopy ( [Mirabel et al. 1991 ; Vanzi et al. 2002 ). 


IRAS 19254-7245 was previously observed in hard X- 
rays with ASGA (Imanishi et al. 1999; Pappa et al. 200d| ). 
The ASGA data were consistent with two possible power- 
law (PL) models: an unabsorbed very flat model (L ~ 1) 
model or an absorbed (A^h ~ 10^^cm“^) and steeper 
(r ~ 1.7) model. No Fe-Ka lines in the range 6.4-7 keV 
were clearly detected, but the upper limit on the EW 
of such lines was ^ 1-2 keV. The statistic of the ASGA 
spectrum prevented those authors from distinguishing be¬ 
tween a Compton-thick (i.e. reflection-dominated) and 
a Compton-thin (i.e. transmission-dominated) source. 
Moreover, given the ASCA’s modest angular resolution 
(~ 2'), it was not possible to pin-point the optical coun¬ 
terpart of the bulk of the X-ray emission. 


3. Data Reduction 


IRAS 19254-7245 was observed by 'XMM-Newton in 
March 200 1 with the EPIC ( Europe an Photon Imaging 
Camera: Btriider et al. 20dl| and Turner et al. 2001 ) 
cameras operating in full-frame mode and the thin filter 
applied. Data have been cleaned and processed using 
the Science Analysis Software (SAS 5.3) and analyzed 
using standard software packages (FTOOLS 5.0, XSPEC 
11.0). The latest calibration files released by the EPIC 
team have been used. Event files produced from the 
pipeline have been filtered from high-background time 
intervals and only events corresponding to pattern 0-12 
for MOS and pattern 0-4 for PN have been used (see the 
XMM-Newton Users’ Handbook Ehle et al. 2001). The 


net exposure time, after data cleaning, are ~ 14.80 ksec, 
~ 18.38 ksec, and ^ 18.39 ksec for PN, MOSl, and MOS2 
respectively. 


IRAS 19254-7245 is detected with a S/N ratio greater 
than ~20 in all the EPIC cameras. In Fig. 1 we show the 
3-10 keV contour plots from the MOS2 image overlaid on 
the DSS2 optical image of IRAS 19254-7245. The hard 
X-ray emission comes mostly from the southern Seyfert 
2 nucleus 0. Since in this paper we want to address the 

^ The soft (E< 2 keV) X-ray emission, too, comes mostly 
from the southern nucleus. The soft X-ray emission appears to 
be extended and it inclndes the whole system. A detailed spec¬ 
tral and imaging analysis of this component will be reported 
elsewhere (Franceschini et al. 2002). 


X-ray spectral properties of the nuclear X-ray emission, 
X-ray counts have been extracted from a circular region 
of 15 arcsec radius^ positionally coincident with the core 
of the hard X-ray emission. Background counts have been 
extracted from a closeby source-free circular region. We 
have then generated our own response matrices (that in¬ 
clude the correction for the effective area) using the SAS 
tasks arfgen and rmfgen. 

In order to improve the statistics, the MOSl and 
MOS2 data have been combined together and the com¬ 
bined MOS and PN spectra have been htted simulta¬ 
neously keeping the relative normalizations free. All the 
models discussed in this paper have been filtered through 
the Galactic absorbing column density along the line of 
sight to IRAS 19254-7245 (A^h = 5.95 x lO^Ocm-^). Unless 
otherwise stated, in what follows errors will be given at 
the 90% confidence level for one interesting parameter 
(Ax2=2.71). 



CENTER: R.A. 19 31 21,44 DEC -72 39 22.1 

Fig. 1. DSS2 image (1' x 1') of the field centered on IRAS 
19254-7245. Contours of the hard (3-10 keV) X-ray emis¬ 
sion have been overlaid on the optical image. We have 
used here the MOS2 data since this detector has the best 


point spread function (see Ehle et al. 2001 ). The displayed 
X-ray contour levels correspond to 5cr, 7a, 10a, 20cr, 30cr 
above the background. 


4. Spectral Analysis 


Single-component models, e.g. a thermal model or 
a PL model, are both rejected by the data at a high 
confidence level. In Fig. 2 we show the ratio between the 
best-fit PL model and the data: large excesses are clearly 
present in the soft energy band (suggesting the presence 


^ The PSF of XMM-Newton does not allow us to extract the 
spectra from a smaller region. 

® The ther mal emissions used in this paper are described by 
the MFKAL (Mewe et al. 1985) model in the XSPFC package. 
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of a soft thermal component) and a “line like” feature 
is evident between 6 and 7 keV (suggesting the pres¬ 
ence of an Fe line). On the other hand a single thermal 
component cannot account for the emission at E > 2 keV. 


We then tried to fit the data with a thermal compo¬ 
nent, an absorbed PL component plus a Gaussian line 
around 6.4 keV. The best-fit values of this model are: kT = 


0.86;° J® keV, r = 1.3+^;^^, A^h = 3.1+^;^ x lO^icm"'' 
(model A in Table 1 and Table 2). The energy position 
and the EW of the Gaussian line are E = 6.48^q 
and EW = keV (source rest-frame). The energy 

position of this line is consistent with cold Fe-Ka emis¬ 
sion. It is worth noticing that the line profile appears to 
be marginally broad, with a width a = 0.31q' 2 keV. This 
broadening can be due to the merging of several lines or 
to a more complex modeling of the continuum (e.g. the 
presence of an absorption edge modifying the underlying 
continuum, see below). 

The flat photon index, F ~ 1.3, measured for 

the hard X-ray continuum, and the presence of a Fe- 
Ka emission are both typical of high-mass X-ray bi¬ 
naries (HMXB) spectra ( White et al. 1983 ). However, 
the equivalent width of the Fe-K complex arising from 
HMXBs, either observed directly in HMXBs (typically 
EW 0.3 keV: see White et al. 1983) or inferred 
for galaxies with X-ray emission dominated by HMXBs 
(Persic & Rephaeli 2002), is not consistent with the huge 
value observed here. The flat hard continuum and the 
prominent Fe-Ka emission line, may alternatively sug¬ 
gest that IRAS 19254-7245 hosts a Compton-thick AGN 
(e.g. Maiolino et al. 19981 ). In this perspective we have 
tried a composite model which accounts for the SB emis¬ 
sion and for the “putative” Compton-thick AGN emission. 
The SB emission has been modeled with a soft thermal 
component. The AGN emission has been described W a 
pure Co mpton-reflected continuum (m odel PEXRAV Q in 
XSPEC, Magdziarz fc Zdziarski 1995 ) combined with an 
absorbed PL model and a Gaussian line at ^ 6.4 keV. 
In this model the absorbed PL component represents the 
scattered emission of the central source, hence its photon 
index is tied to the one in the PEXRAV model. The re¬ 
sults are reported in Tables 1 and 2 (model B) and shown 
in Fig. 3. We note that the best-fit photon index is ~1.8, 
a value very similar to that observed in unobscured AGN 
( [George et al. 2000 ). 

As described above, the baseline model for the hard X- 
ray continuum of IRAS 19254-7245 comprises a reflected 
plus a scattered component. If this latter component is 
supposed to arise from a “warm mirror”, then one would 
also expect to observe emission lines from highly ionized 
Fe. We have estimated the upper limits on the EW of 


Since we lack any information about the emission at energy 
above 10 keV we kept the inclination, the abundances and the 
energy cutoff frozen to standard values of ~ 60°, 1 and 100 
keV; only the intrinsic PL photon index and the reflection nor¬ 
malization were allowed to vary in the fitting procedure. 


the He-like (6.7 keV) and the H-like (6.97 keV) Fe emis¬ 
sion lines to be ~ 1.8 keV and ~2 keV, respectively. A 
second narrow gaussian line could probably be accommo¬ 
dated within the data uncertainties. These upper limits 
are consistent with that expected from the scattered com¬ 
ponent, but the present statistics does not allow us to be 
conclusive about this point. 

In summary, the broad-band X-ray emission from 
IRAS 19254-7245 requires a thermal emission {kT ^ 0.85 
keV) plus a hard continuum. The spectral parameters of 
the thermal emission are not a strong function of the de¬ 
tailed modeling of the hard X-ray continuum; this com¬ 
ponent (which dominates at A ^ 1 keV) is likely associ¬ 
ated with the starburst and, compared to the hard one, 
appears to be spatially extended. A detailed analysis of 
this emission will be presented in a forthcoming paper 
(Franceschini et al. 2002). The hard X-ray continuum can 
be described by two models: a) a Compton-thin (Table 1 
model A) AGN having a very flat photon index or b) a 
Compton-thick (Table 1 model B) AGN having an intrin¬ 
sic photon index of ~ 1.8. Both models require a strong Fe 
line at ~6.49 keV (rest-frame). Since these two competing 
models give similar ')^ jv in the next section we will try to 
discriminate between them. 


data/model 


1 "M" J III 1 tl-4-^ 

Jj-f “F 


' 1 ' 


. .1_1_ ■ 

0.5 1 2 5 

channel energy (keV) 


Fig. 2. Ratio between the PN (red) and MOS (black) data 
and the model when the spectrum of IRAS 19254-7245 is 
fitted with a single power law model; strong residuals are 
present at energy below 1 keV and above 5 keV. 


5. Discussion 

A viable interpretation of the observed flat (F ~ 1.3) 
hard X-ray continuum and prominent {EW ^1.4 keV) 
Fe-Ka emission line involves a reprocessing of the pri¬ 
mary spectrum, i.e. the presence of a Compton-thick 
source. Further indications of the Compton-thick nature of 
IRAS 19254-7245 come from observations at other wave¬ 
lengths. One of the broad-band diagnostics largely used 
in the literature to estimate the strength of the nuclear 
source and the Compton thickness is the ratio between 
the strength of the [ 0777 ]A 5007 emission line and the 2 - 
10 keV flux. This ratio -once that the [O 777 ] emission 
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Table 1. Results from the best-fit model. 


Model 

KT 

Norm“ 

r 

Ah 

Norm" 

Norm" 


EW 

xVdof 


keV 






keV 



A 

B 

0.86:";i^9 

0.85t°;J^ 

1 Qqi-U.CJU 
i.zo_0.3i 

i.^u o.39 

1.0U_Q 21 

1 

0 Q1 

U.Ol_0.i5 

0 47+0-52 
-0.25 

2-63:S"I 

q y-i +2.40 
‘ -*--1.40 

/ 

6.48:";- 

6.491^® 

1.4±0.5 

2.0±0.6 

50.9/57 

48.1/56 


“ in units of 10“^®/(47r(D^) f HeTindV , where D is the distance to the source (cm), Ue is the electron density (cm“^), and nn 
is the hydrogen density (cm“®) at 1 keV. 

** in units of 10“® photons keV“^ cm~^ at 1 keV; in model B it refers to the scattered component, 
in units of 10“® photons keV“^ cm“^ at 1 keV; in model B it refers to the reflected component. 

Table 2. Observed X-ray Fluxes and X-ray luminosities 


Model 

BAND 

10-14 

Flux 

erg cm“^ s“ 

•1 

Luminosity" 

10'‘^ erg 



TOTAL 

Starburst 

AGN 

Starburst 

AGN 

A 

0.5-2 keV 

5.30 

2.13 

3.16 

4.4“ 

9.7" 


2-10 keV 

23.0 

0.1 

22.9 

0.30“ 

37.7" 

B 

0.5-2 keV 

5.31 

2.10 

3.21 

4.3“ 

13.6" 


2-10 keV 

23.5 

0.1 

23.4 

0.28“ 

38.7" 


“ X-ray luminosity de-absorbed from intrinsic and galactic Ah. 
^ X-ray luminosity de-absorbed from galactic Ah. 



0.5 1 2 5 10 

channel energy (keV) 


Fig. 3. IRAS 19254-7245 (Model B). Upper panel: un¬ 
folded X-ray spectrum. Lower panel: ratio between the 
data (MOS-I-PN) and the best-fit model. 


is corrected for the extinction using the Balmer decre¬ 
ment ( Bassani et al. IQ'Q^ - combined with the EW of the 
measured Fe emission line can be used to evaluate the 
Compton thickness of the source. 

Using recent [O///] estimates (for the southern nu¬ 
cleus; see Berta et al. 2002) and the 2-10 keV observed flux 
(see Table 2) we obtain F(2-10 keV)/F[ 0777 ] ~ 0.09. This 
value, combined with the measured EWpe-Ka ~ 1 — 2 
keV of the Fe line, clearly locates IRAS 19254-7245 in the 
region populated by Compton-thick sources in Fig. 1 of 


Bassani et al. 1999 . 

The energy position of the line (6.49tQ °o keV) is con¬ 
sistent with cold (i.e. Fe less ionized than Fe XVII) Fe Ka 
emission but not with highly ionized (He-like or H-like) 
iron. Therefore we can exclude an origin of the line from 


optically thin gas (the so called “warm mirror”) placed 
outside the torus (see e.g. Matt et al. 1996| ), leaving the 
production from a cold medium as the only possibility. In 
this context the line could be produced by reflection from 
the inner surface of the circumnuclear absorbing mate¬ 
rial ( IChisellini et al. 1994) or by trans mission through the 
torus itself ( Leahy fc Creighton 19931 ). Moreover the X- 
ray spectrum is well described by a reflection component, 
and the equivalent width of the detected line (with re¬ 
spect to the reflected component) is 2 keV: this is in good 
agreement with the expectation from a reflection domi¬ 
nated source (see Matt et al. 1996 ). 

If the Compton-thick hypothesis is correct, then the 
observed luminosity (see Table 2) could be just few 
percent of the intrinsic one. Since the best-fit luminos¬ 
ity of the scattered component is 

then the intrinsic luminosity is estimated to be 
= Assuming the NGC1068 


L 


INT 


(2-10) 


values ( Iwasawa et al. 1997 ) for the mirror optical depth 
(r ~ 10“^) and the subtended angle (n/27r = 0.25), we 
obtain an intrinsic luminosity for the AGN harbored in 
IRAS 19254-7245 of ~ lO^^erg s"!. 


The derived intrinsic luminosity of the AGN present 
in IRAS 19254-7245 suggests that the AGN activity 
contributes a significant fraction of the FIR luminosity. 
Indeed mid-IR observations suggest that the AGN con¬ 
tribution to the FIR emission could be 40% (see Berta 
et al. 2002). We have then computed log( Lfir/A 2 -io) 
for the AGN component. We found that this ratio is 
^ 3, in agreement with the Compton-thick scenario 
(Mulchaey et al. 1994). 

In Fig. 4 we compare the spectral energy distribution 
(SED) of IRAS 19254-7245 with the SED of NGC6240 
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(Iwasawa et al. 2001). The latter ULIRG is known to har¬ 
bor a heavily-obscured, high-luminosity AGN. In this re¬ 
spect it is worth emphasizing the similarities of the two 
SEDs from the IR-to-X-ray band: in particular, both 
sources show a strong Fe emission line superimposed on a 
strongly absorbed continuum. In more detail NGG6240 
and IRAS 19254-7245 show a similar emission in the 
X-ray band, while they differ in the optical/FIR band 
(IRAS 19254-7245 being more luminous then NGG6240). 
This could indicate that the AGN in IRAS 19254-7245 
could be intrinsically more powerful and/or more ob¬ 
scured in the X-ray energy band. Alternatively, the SB 
activity may be stronger in IRAS 19254-7245 than in 
NGG6240. Finally the shapes of these SEDs are simi¬ 
lar to the SED of the cosmic energy density spectrum 
(see Lehmann et al. 2001 and reference therein). In both 
ULIRGs there is evidence of comparable contributions by 
obscured SB and AGN activities. 



Fig. 4. Comparison between the SED of the Southern nu¬ 
cleus of IRAS 19254-7245 (open squares) and NGC6240 
(filled squares and solid line). The data of the SED of 
NGC 6240 are from llwasawa et al. 20011 The data of 


IRAS 19254-7245 are from Berta et al. 2002. 


6. Conclusions 

We have presented the XMM-Newton observations of 
IRAS 19254-7245. The main results are hereby summa¬ 
rized. 

1 . The hard X-ray emission originates from the southern 
nucleus of IRAS 19254-7245, which is optically classi¬ 
fied as a Seyfert 2 galaxy. 

2. We detect, for the first time in this system, a strong 
{EW ~ 1.4 — 2.0 keV) emission line at 6.49 ± 0.10 keV 


(rest-frame), which is consistent with a cold Fe-Ka 
line. 

3. The broad-band X-ray spectra requires a soft thermal 
component, which accounts for the starburst emission, 
and a flat (T ~ 1.3) continuum most likely associated 
with a heavily absorbed AGN. The hard continuum 
and the detected strong Fe-Ka emission line are highly 
indicative of a Compton-thick source. Indeed, from a 
statistical and a physical point of view, the X-ray spec¬ 
tral properties of this object can be well explained with 
a Compton thick model. In this scenario, the intrin¬ 
sic 2-10 keV luminosity of IRAS 19254-7245 could be 
greater than 10"^'^ erg s“^, i.e. attaining the QSO lu¬ 
minosity regime. 
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